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Abstract 

A numerical model for the solar modulation of cosmic rays, based on the solution of a 
set of stochastic differential equations, is used to illustrate the effects of modifying the 
heliospheric magnetic held, particularly in the polar regions of the heliosphere. To this 
end, the differences in the modulation brought about by each of three choices for the 
heliospheric magnetic held, i.e. the unmodihed Parker held, the Smith-Bieber modihed 
held, and the Jokipii-Kota modihed held, are studied. It is illustrated that both the 
Jokipii-Kota and Smith-Bieber modihcations are ehective in modifying the Parker held 
in the polar regions. In addition, it is argued that the modihcation of Smith and Bieber 
is based on observational evidence and has a hrm physical basis, while these motivations 
are lacking in the case of the Jokipii-Kota modihcation. From a cosmic ray modulation 
point of view, we found the Smith-Bieber modihcation to be the most suitable choice for 
modifying the heliospheric magnetic held. The features and ehects of these three mod¬ 
ihcations are illustrated both qualitatively and quantitatively. It is also shown how the 
Smith-Bieber modihed held can be applied in cosmic ray modulation models to reproduce 
observational cosmic ray proton spectra from the PAMELA mission during the solar mini¬ 
mum of 2006 - 2009. These results are compared with those obtained in previous studies of 
this unusual solar minimum activity period and found to be in good qualitative agreement. 

Keywords Solar modulation ■ Cosmic rays ■ Stochastic diherential equations ■ Helio¬ 
spheric magnetic held ■ PAMELA ■ Solar minimum 

1 Introduction 

The intensities of Galactic cosmic rays (CRs) entering the heliosphere are modulated by a 
number of physical mechanisms, including convection, particle drifts and dihusion, as well 
as adiabatic energy changes. Modeling these mechanisms and the associated modulation 
of CRs is done primarily by employing numerical modulation models, since analytical 
models are limited to only a few simplihed cases. 

The numerical modulation model used in this study is based on the solution of an 
appropriate set of stochastic diherential equations (SDEs). Studies in the past have mainly 
employed numerical models based on hnite diherence approaches, the alternating direction 
implicit (ADI) being the most extensively utilised scheme. However, after Zhang (1999) 
successfully applied SDEs to study CR modulation in the heliosphere, the use of such 
models has become popular. SDE-based models have been used to study the transport 
of pickup ions (Fichtner et ai, 1996), solar energetic particles (Droge et ai, 2010), the 
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propagation of Galactic and Jovian electrons (Strauss et ai, 2011a), and other aspects 
of CR modulation (e.g. Pei et al. 2010; Strauss et ai, 2011b, 2012; Luo et ai, 2011, 
2013a, 2013b). Some historical and contemporary models are discussed by e.g. Yamada, 
Yanagita and Yoshida (1998), Alanko-Huotari et al. (2007), Bobik et al. (2012), and 
Kopp et al. (2012) who also discussed the mathematical and practical implementation of 
SDEs for a variety of general transport equations. 

The SDE-approach leads to numerical models that are ideal for implementation on 
multiple processors, allowing for parts of the numerical code to be executed simultane¬ 
ously and independently, saving a great deal of computation time. Besides the parallel use 
of multiple central processing units (CPUs) and computer clusters, high performance cal¬ 
culations on graphics processing units (GPUs) have also become available during the last 
years; in such a GPU-accelerated implementation to study the solar modulation of CRs, 
presented by Dunzlaff, Strauss and Potgieter (2015), a performance increase of a factor of 
about 10 - 60 was found when compared to the CPU-version of the same algorithm. ADI- 
based models, on the other hand, only allow for a linear execution of the numerical code so 
that improvement in computation time cannot be achieved by these means. Furthermore, 
SDE-based models, being independent of a spatial grid, have also displayed remarkable 
numerical stability. This makes it possible to incorporate detailed descriptions of the 
heliosphere and heliospheric structures into these models. Examples of such implemen¬ 
tations include the works of Strauss et al. (2011b), where the relative motion of Jupiter 
— as the source of Jovian electrons — was accounted for; Strauss et al. (2012), where 
a fully three-dimensional wavy heliospheric current sheet (HGS) was implemented; Luo 
et al. (2013a), where acceleration of the anomalous component and the re-acceleration 
of Galactic GRs by the solar wind termination shock (TS) were included; and Luo et al. 
(2011, 2013b), where time-varying modulation parameters were implemented. ADI-based 
models, on the other hand, are notoriously unstable when solving differential equations 
in higher dimensions and the complexity of these models are thus restricted. In addition 
to these considerations, SDE-based models are especially powerful in their ability to vi¬ 
sualise the modulation process, being able to calculate pseudo-particle trajectories or — 
equivalently stated — follow the evolution of individual phase space (density) elements; 
it is important to note that these trajectories are not the trajectories of actual particles, 
but can at most be interpreted to obtain an indication of how the modulation of actual 
particles would proceed. The ability to calculate these trajectories leads to the possibility 
of calculating additional modulation features which were previously not possible, e.g. the 
propagation times and energy losses of particles in the heliosphere. 

The heliospheric magnetic held (HMF) plays a central role in the modulation mecha¬ 
nisms pointed out above and, as such, the choice of the HMF prohle is a very important 
question in modulation studies. In this paper, we apply the SDE-based model to investi¬ 
gate the effect of different choices for the HMF. We consider the unmodihed Parker HMF 
(Parker, 1958) and two modihcations to this held, namely that of Jokipii and Kota (1989) 
and that of Smith and Bieber (1991). The focus is on the modihcation of Smith and 
Bieber (SBM), and it is argued that this modihcation should be preferred when modify¬ 
ing Parker like HMFs. To this end the numerical model is applied, employing the SBM, to 
reproduce the Galactic proton spectra observed at Earth by the PAMELA mission during 
the peculiar solar minimum of 2006 to 2009. The results are then compared to those of 
other authors who have studied this minimum. 
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2 The transport model for Galactic protons 


An equation including most of the relevant heliospheric modulation mechanisms was de¬ 
rived by Parker (1965) and is known as the Parker transport equation (TPE), given by 

^ = -(vS„ + {fl'd>)-v/ + v.(K..v/) + t(v.C„)^ + (3, (1) 


in terms of the omnidirectional distribution function / {f,p,t), where Pis the position, p 
the particle momentum, and t the time. The hrst term on the right hand side describes 
the process of convection via the expanding solar wind (SW), with velocity the second 
term represents the effects of gradient and curvature drifts, where {v^) is the pitch-angle 
averaged guiding centre drift velocity; the third term describes particle diffusion through 
the symmetric diffusion tensor Kg; the fourth term accounts for adiabatic energy changes; 
the hnal term on the right hand side, Q, includes any source terms, if required. For a 
detailed description of these modulation processes, see e.g. Potgieter (2013). 

The geometry, magnitude and turbulence of the HMF play a key role in most of these 
modulation processes, in particular in determining particle drifts, as will be shown and 
emphasised below. In most CR modulation models the straight-forward Parker HMF 
(Parker, 1958) is used. This HMF prohle is valid for r > r©, where r© = 0.005 AU is the 
solar radius, and is expressed in spherical coordinates (r, 9, 0) as 
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where and are unit vectors in the radial and azimuthal directions respectively, Bq is 
the magnitude of the HMF at tq = 1 AU (i.e. at Farth), and the spiral angle "0 is dehned 
by 

= = (3) 

^SW 

In the above expression, H = 2.66 x 10“® rad.s“^ is the average angular rotation speed of 
the Sun. The magnitude of the Parker HMF is then given by 
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Figure shows the magnetic held lines in the case of the Parker HMF, converted to 
Cartesian coordinates, and its spiral structure is clearly illustrated over the hrst 10 AU 
from the Sun, which is located at (X, Y, Z) = (0, 0, 0); these lines originate at consecutive 
azimuthal angles and a latitude oi 9 = 45°. 

The Parker HMF, however, leads to very large gradients in particle intensities j = fP^, 
with P the particle rigidity, in the polar regions; this makes for a gross over-estimation 
of drift ehects in these regions, as was originally pointed out by Potgieter, le Roux and 
Burger (1989); see also Potgieter (2013). An example of such drift dominated modulation 
is found in Jokipii and Kopriva (1979). The Parker HMF has therefore been modihed 
in several studies to scale down these ehects in the polar regions; this forms the topic of 
Section HI 

The SW prohle selected is rather simple, accelerating from zero to a constant and 
supersonic velocity within the hrst 0.3 AU from the Sun. The only other radial structure 
occurs at the position of the TS, i.e. at r = tts, where Uw drops oh to subsonic levels; 
The SW also has a latitudinal prohle, which accounts for a slow SW of ~400 km.s“^ in 
the equatorial regions and a fast SW of ~800 km.s“^ in the polar regions, in accordance 
with Ulysses observations during solar minimum (Phillips et ai, 1995). 
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Particle diffusion is realised through the diffusion tensor Kg in Eq. Q, of which each 
element is known as a diffusion coefficient. An approach similar to that of Potgieter et 
al. (2014) is followed; this will also allow for a direct comparison with their modulation 
modeling, which did not use the SBM. The parallel diffusion coefficient k\\, i.e. parallel 
to the mean HMF, is assumed to have a power-law dependence on particle rigidity, while 
its spatial dependence is assumed to be inversely proportional to the magnitude of the 
magnetic held, so that 


0,2 — 


where is a constant in units of 6 x 10^° cm^.s“^, with Pq = 1 GV and Bq = 1 nT 
added to obtain the correct dimensions; (3 = v/c is the ratio of particle speed v to the 
speed of light c. Here Ui and 02 are dimensionless constants and determine the slope of 
the rigidity dependence below and above a rigidity Pk respectively; in this study, we set 
02 = 1.95. The quantity 03 = 3.0 is another dimensionless constant and determines the 
smoothness of the transition between the two slopes P“^ and P“^ at Pk. The rigidity 
dependence is therefore essentially a double power-law. 

The perpendicular diffusion coefficient kj_ consists of perpendicular diffusion in the 
and 60 directions respectively, i.e. K±r and k± 0 . Here, k,± is scaled as «:||. Giacalone and 
Jokipii (1999) found that the ratio n±/n\\ has a value between 0.02 and 0.04. In addition, 
observations from the Ulysses spacecraft have revealed that the latitude dependence of GR 
protons is signihcantly less than predicted by classical drift models (Potgieter and Haas- 
broek, 1993), which led Kota and Jokipii (1995) to propose the concept of an anisotropic 
perpendicular diffusion, where in the off-equatorial regions (e.g. Burger, Pot¬ 

gieter and Heber, 2000; Ferreira et al, 2000; Potgieter, 2000). This anisotropy of is 
accounted for in this work, as was done by e.g. Ngobeni and Potgieter (2008) as well as 
Potgieter et al. (2014), so that 
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where = 0.02 are dimensionless constants, and 
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where d = 3.0 is a dimensionless constant that determines the enhancement factor of 
K ,±0 from its value in the equatorial plane towards the poles, with respect to kh; see also 
Potgieter (2000). 

Cosmic rays will undergo a combination of gradient and curvature drifts (ud)gc caused 
by the large scale HMF, as well as current sheet drift (hd)^^ because of a switch in magnetic 
polarity over the HCS. These drifts are both calculated from the general equation for the 
pitch-angle averaged guiding center drift velocity, written as 
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and = B/B; H is the Heaviside function, ^Dirac is the Dirac delta function, 9' is the 
HCS latitudinal extent dehned below, and ka the global drift coefficient which is given 
by 




PP jP/Pof 
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(13) 


In this equation q is the particle charge sign so that in this study, which only considers 
results pertaining to CR proton modulation, q = -|-1. The quantity Pq G {1/vTo, 1/V^} 
GV is added on dimensional grounds. The value A = ±1, i.e. either H > 0 or H < 0; this 
represents the magnetic polarity cycle with H > 0 (H < 0) dehned as the polarity cycle in 
which the HMF points outward (inward) in the Northern hemisphere of the Sun; during 
an H > 0 (H < 0) polarity cycle positively charged particles will drift inward (outward) 
along the polar regions and outward (inward) along the HCS in the equatorial regions, 
while the reverse is true for negatively charged particles. The quantity Ka is a constant 
which can be used as a scaling factor, having values between 0.0 (zero drift) and 1.0 (full 
drift); in this work, = 1.0, which is known as the so-called weak scattering maximal 
value for the drift coefficient. Remembering the change in sign that occurs when crossing 
the HCS, the gradient and curvature drifts can thus be calculated by means of 
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where Kp^ was dehned as 


Ka = Ka^H 
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(15) 
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For the Parker HMF, the averaged gradient and curvature drift velocities in the radial, 
latitudinal, and azimuthal directions respectively are given by 


= —Udo cot 9 
= Udo(2 + d'2) 
= Udod' cot 6*, 
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Figure shows the drift velocity streamlines for this case, i.e. solving for 
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where f = (r, 0, 0) is the position vector in spherical coordinates and I is the arc length 
along the held lines of (hd)gc. These lines are calculated in the polar regions and it is 
clear from their narrow, cone-like conhguration that drift in these regions occurs very 
effectively, towards the Sun. Note that Figure again makes use of Cartesian coordinates 
so that the Sun is located at {X,Y,Z) = (0,0,0); the polar regions are centered around 
the Z-axis. 

When accounting for current sheet drift, which is directed parallel to the HCS and 
perpendicular to the HMF (e.g. Burger, Moraal and Webb, 1985; Burger and Potgieter, 
1989), the model does not make use of the ^oeita function in Eq. 0 ; rather, it is assumed 
that a particle will experience current sheet drift upon the condition that its distance d^s 
to the HCS is less than two gyro radii, i.e. 


dns < 2rL = 2 
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with q the particle charge (not to be confused with the particle charge sign q as in Eq. 
(13)), and m is the relativistic mass of the particle. Assuming that the latitudinal extent 
of the HCS is given by the expression from Kota and Jokipii (1983) 
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Eq. (16) is replaced by 


sin -0 cos ( ±/3rot) 

Vde = Vns,0 sin (±/3rot) 

Vd^ = Vns,oCOS^p cos ( ±f3,ot ) , 


( 21 ) 


when the condition of Eq. (19) is met. Here Uns,o is given by the approximation of Burger, 
Moraal and Webb (1985) 
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and /drot is the angle between the radial direction and a vector parallel to (Fd)ns> dehned 
by 

do' 

tan/3rot = r —, (23) 

as in Strauss et al. (2012) and discussed by Burger (2012). The correct sign of /3 in Eq. 
(21) is determined by 

do' 

— <0 ^ /?rot < 0 

dr 
dO' 

— >0 ^ Prot > 0 . 

dr 


( 24 ) 
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The proton LIS assumed in this study has its origin in that of Langner and Pot- 
gieter (2004). Potgieter et al. (2014) improved on this LIS by taking into account 
PAMELA measurements at energies between 30 and 50 GeV. However, this LIS did not 
provide for the heliopause (HP) crossing of Voyager 1 and was therefore modihed to take 
into account these new Voyager 1 observations in the lower energy range; see Potgieter 
(2014) in this regard. The LIS resulting from these modihcations, expressed in units of 
particles.m“^.s“^.sr“^.MeV“^, is employed in this study and is given by the expressions 

j = 0.6978 exp {4.64 - 0.023(ln Ef - 2.91E-0-5} 
for energies E < 1.4 GeV, and 

j = 0.6847exp {3.22 - 2.78(lnE) - 1.5^"^} (25) 

for E > 1.4 GeV. The HP position is assumed at 120 AU, where this LIS is specihed as 
an initial condition for cosmic ray modulation. 


3 The Numerical Modulation Model 


Since analytical solutions to the TPE in Eq. ([^ are only possible to a very limited extent, 
modulation models solve this equation numerically to various degrees of complexity. The 
numerical model of this work makes use of an SDE approach and is based on the well- 
benchmarked model of Strauss et al. (2011a, 2011b, 2012). It is only briefly explained 
here; an extensive treatment can be found in Kopp et al. (2012). The model is solved in a 
steady state and assuming no sources so that Q = 0. It is solved in spherical coordinates 
(r, 6, (p) and particle kinetic energy E. 

Each Fokker-Planck type equation has a corresponding set of independent SDEs, which 
is calculated in either a time forward or time backward fashion. The TPE, written in 
spherical coordinates, is associated with such a time backward set of SDEs, given by 
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(26) 

for each of the spherical coordinates r, 9, and 0, as well as for the particle energy E; 
dW = [dWr, dWe, dW^] is the multi-dimensional Wiener process with each of its elements 
containing a Gaussian distributed random number. The quantity P is defined by 


E -|- 2Eq 
E + Eq 


( 27 ) 


where Eq is the particle rest energy; ds is the inhnitesimal backward time increment. 
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The diffusion tensor has here been converted to spherical coordinates and also in¬ 
cludes the drift coefficient ka. Its elements are related to the coefficients in HMF aligned 
coordinates through 
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(28) 

This tensor transformation is only valid for Parker-like HMF profiles, i.e. for an HMF 
that does not contain a 6 component; more general transformations are discussed in e.g. 
Effenberger et al. (2012). 

Solving the TPE in a time backward fashion implies that a number of pseudo-particles, 
iV, are traced from Earth outward to the boundary of the heliosphere, stepping in (r, 6^, </>) 
and calculating energy losses according to Eq. (26). At the outer modulation boundary, 
?"out = 120 AU, particle relative intensities are convolved with the LIS and an energy 
spectrum at Earth is obtained. For the purposes of this study, no modulation is assumed 
beyond the HP, so that the heliopause spectrum and the very LIS is synonymous (see 
Potgieter et al., 2013). A reflective inner modulation boundary is used at rjn = 0.1 AU, so 
that any pseudo-particle for which r < rjn is reflected back into the modulation volume, 
i.e. Hn < r < rout- 


4 Modified Magnetic Field Profiles 


4.1 The modification of Jokipii and Kota 


Jokipii and Kota (1989) showed that the effect of small perpendicular magnetic held com¬ 
ponents near the solar surface leads to a larger magnetic held at greater radial distances 
in the polar regions. These authors argued that the average direction of these small com¬ 
ponents would cancel out so that only the magnitude of the magnetic held is modihed. 
They put forward the following modihed expression for the Parker HMF magnitude 


B 




(29) 


introducing the quantity 6 to signify the magnitude of the extra, superimposed, transverse 
magnetic held. The Jokipii-Kota modihcation is henceforth referred to as the JKM, while 
the term containing 6 will be referred to as the JKM modihcation term. In this study we 
use 6 = 0.002 (see e.g. Potgieter, 2000). Geometrically speaking, a cone centered around 
the polar regions is induced in which the HMF is altered so that it decreases as r~^ instead 
of r“^. The ehect of this modihcation, therefore, is to bring about the changes which are 
required to reduce the large drift velocities in the polar regions, without noticeably altering 
the held in the equatorial plane. 

However, upon this modihcation, the requirement that the magnetic held remains 
divergence free, i.e. V ■ B = 0, is no longer complied with and therefore the JKM 
is technically incorrect. Through inspection, the HMF will remain divergence free if 
J = 6{9) (X (sin6*)”^, so that 


S(B) = 


(30) 


sin^’ 

with Jm = 8.7 X 10“® (Langner, 2004). This value for Jm implies 6{6) = 0.002 near the 
poles and J(0) ~ 0 in the ecliptic plane. This will again bring about the intended changes 
to the HMF. 











4.2 The modification of Smith and Bieber 


This modification was introduced by Smith and Bieber (1991), who analysed data col¬ 
lected by various satellites and found magnetic field spirals at Earth more tightly wound 
than that predicted by the Parker theory. Consequently, they suggested a modification 
to the Parker spiral by reasoning that the difference in rotational speed of the equatorial 
and polar regions of the Sun would cause small azimuthal magnetic field components to 
develop. This would lead to larger spiral angles at larger radial distances, and would 
provide an explanation for the larger spiral angles at Earth. At the same time, this modi¬ 
fication would have a significant effect on the HMF structure at the high latitude regions. 
The expression for the spiral angle, i.e. Eq. ([^, was modified to yield 


il{r — b)sm9 rVsw{b,0) /BT{b)\ 

Kw(r,0) bV,^{r,e) UrW; ’ 


(31) 


where B^ib)/B-£i{b) is the ratio of the azimuthal to the radial magnetic field components 
at a position b near the solar surface, here taken as b = 20rQ. Smith and Bieber (1991) 
showed that the value of B-^ib)/B'^ib) is approximately -0.02 and although the value of b 
assumed here is larger than the original b = 5rQ assumed by these authors, the value of 
-0.02 is retained for the purposes of this paper, except where explicitly indicated other¬ 
wise. This modification, based on sound physical arguments, has been mostly ignored in 
the numerical modeling of the solar modulation of cosmic rays, which is rather surpris¬ 
ing. It will be shown next that this modification is indeed appropriate and as such an 
improvement. 

Figure is an illustration of HMF lines that results from applying the SBM; these 
lines once again originate at consecutive azimuthal angles and at a latitude of 6* = 45°. 
According to Eq. (31), and clearly evident from comparing Figures and the SBM 
alters the direction of the Parker HMF by winding up the magnetic field spirals more 
tightly; as a consequence, in contrast to the straight line case of Eqs. ([^ and (Q, the 
HMF still has an azimuthal component in the polar regions, even at 6 = 0°. It is also 
clear from this comparison that the field lines in the case of the SBM are closer together 
than in the case of the Parker HMF, in accordance with the larger HMF magnitudes of 
the modified field. Note that, in Figure]^ B-^ib)/B-p^ib) = —0.5 was chosen to be very 
small, corresponding to a very large modification, in order to illustrate these effects more 
clearly. 

It should be noted that the ultimate, very complex, HMF modification was introduced 
by Fisk (1996). However, it has been seen as quite controversial and without convincing 
observational support so that it is not pursued any further in this work; see Sternal et 
al. (2011) for an appreciation of this modification. The SDE-approach does however 
allow for the implementation of a HMF of this complexity whereas traditional ADI based 
numerical schemes are quite unsuitable. 


4.3 Comparison of the features of HMF modifications 

The top panel of Figure compares the radial profiles of the Parker HMF and HMFs 
which have been modified according to respectively the SBM and JKM; this is shown 
in the equatorial plane with 6 = 90°, assuming a magnetic field magnitude of 5.05 nT 
at Earth. All three of these HMF profiles start out with an dependence at the 

smallest radial distances and end with an dependence just in front of the TS at 

f' = ''"TS) where 14^ drops by a factor of 2.5 and brings about a step-like increase; beyond 
the TS, the dependent decrease is continued. In the case of the Parker HMF, the 
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j,-i.95 (dependence continues up to a radial distance of about 0.5 AU, where it begins to 
transition into a weaker dependence to eventually reach at ~ 3.0 AU. Although not 

clearly discernible on the scale of this graph, the SBM and JKM prohles were determined 
to enter the dependence by a distance of up to 1 AU earlier than the Parker HMF. 

The bottom panel of Figure compares the radial profiles of the Parker HMF, SBM, 
and JKM in the polar regions and necessarily illustrates the effects of the modifications 
more clearly. The Parker HMF and both the SBM and the JKM start out with an 
dependence at the smallest radial distances. At a radial distance between 6 AU and 7 
AU the Parker HMF starts transitioning into a weaker radial dependence and keeps up 
this transition to reach its weakest dependence of at a radial distance of between 40 
AU and 50 AU; this is then maintained up to the TS. The SBM begins its transition into 
the weaker radial dependence at a smaller radial distance between 1 AU and 2 AU and 
reaches at about 20 AU, continuing up to the TS. The JKM begins the transition 
into the weaker radial dependence even earlier than the SBM at ~ 1 AU and reaches 
r-i oo between 8 and 9 AU, once again continuing up to the TS. In all three cases, the 
dependence in front of the TS is regained beyond the TS. 

Notice that, in the polar regions, the JKM profile does not undergo the upward step 
at the TS to the extent shown by the Parker HMF and SBM; in fact, the JKM profile 
barely shows any visible increase. This suppression of the TS in the polar regions is not 
an intended effect and is a consequence of the fact that the TS is simulated entirely by 
the drop in Uw at r = rxs- In the instance of the JKM, the modification term in Eq. 
(29) is much larger than the term containing Uw and hence the effect of the drop in Uw 
at r = rxs is markedly less pronounced. 

It is clear from Figure that both the SBM and JKM have the desired effect of 
drastically reducing the radial dependence of the HMF in the polar regions: the transition 
into the weaker radial dependence is initiated at much smaller r than in the case of the 
Parker HMF, ensuring significantly larger magnetic field magnitudes B. It is also noted 
that because the SBM starts its transition into the weaker radial dependence at larger 
radial distances than does the JKM, its effect on B is less pronounced than that of the 
JKM. Although not illustrated in Figure it is of importance to note that this less 
effectual status of the SBM when compared to the JKM is only locally true, i.e. it is only 
valid at or very close to the poles as in this case with 6 = 5°. 

The four panels of Figure]^ show the unfolding of the Parker HMF, SBM, and JKM 
latitudinal profiles over successive radial distances of 1 AU, 10 AU, 50 AU, and 100 AU. 
From the first panel it is clear that, at 1 AU, the effects of both the SBM and JKM 
are negligible in the equatorial regions, but become more pronounced towards the poles. 
Taking the equator at 6 = 90° as reference and moving towards the poles, the deviation of 
the SBM from the Parker HMF is seen to start at about 6 = 90° ±22°, continuing towards 
the respective poles. The SBM is larger than the Parker HMF up to about 6 = 90° ±76°, 
where it drops below and eventually reaches a maximum deviation from the Parker HMF 
of ~ 0.08 nT at the poles. The difference between the Parker HMF and JKM is seen to 
consistently increase towards the polar regions 6 G {0°, 180°}, where the JKM eventually 
reaches levels of ~ 0.12 nT higher than the Parker HMF. 

The picture at a radial distance of 10 AU has unfolded to be signihcantly different 
from the profiles at 1 AU. The SBM is seen to have already modified the Parker HMF 
in the equatorial regions, amounting to a value of about 0.04 nT higher than the Parker 
HMF at 6^ = 90°. Towards the poles the difference between the SBM and Parker HMF 
increases to about 0.06 nT at intermediate latitudes, declining again to about 0.03 nT 
at 6* G {0°, 180°}. The SBM, at the equator, is initially higher than the JKM, but drops 
below it at about 6 = 90° ± 42°, and continues to do so up to the polar regions. Smaller 
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deviations from the Parker HMF in the eqnatorial regions are also seen to be indnced by 
the JKM, amonnting to no more than 0.02 nT at 6 = 90°. As was the case at 1 AU, the 
effect of the JKM increases towards the poles where it reaches levels of abont 0.10 nT 
higher than the Parker HMF. The pictnre at 50 AU looks qnalitatively the same as that 
at 10 AU, only the magnitnde of the HMF being smaller. Comparing with Fignre it is 
clear that this qnalitative pictnre is continned up to the TS. 

The fourth panel of Figure shows the latitudinal prohles at 100 AU, a distance 
comfortably beyond the TS. The SBM is still higher than the Parker HMF at all latitudes, 
now being higher than the Parker HMF by about 0.01 nT in the equatorial regions and 
by about 0.02 nT in the polar regions. The JKM is, curiously enough, at more or less the 
same level as the SBM in the polar regions; it lies just below the Parker HMF at 6* = 90°, 
but climbs above it at about 6 = 90° ± 45°, eventually reaching the level of the SBM at 
the poles. This unfolding of the HMF latitudinal prohles over successive radial distances 
effectively illustrates the previous statement about the more effective status of the JKM 
being only locally true: the SBM is effected over a larger part of the modulation volume 
so that its overall effect can indeed be greater than that of the JKM. 

From Eqs.(|^ to 0 . dehning the diffusion coefficients straight-forwardly as used in this 
study, it is clear that the choice of HMF prohle directly inhuences the spatial dependence 
of these coefficients. Figure shows the radial prohles of the parallel mean free paths A|| 
of 1 GV protons, which is related to the parallel dihusion coefficient via 


-'ll = “''II' 


(32) 


This is shown for each of the Parker HMF (solid red lines), the SBM (dashed green lines), 
and the JKM (dashed blue lines). The upper three lines represent the values of A|| in 
the polar {6 = 5°) regions, while the lower three, almost coinciding, lines do so in the 
equatorial plane at 0 = 90°. As expected, the ehects of the modihcations are minimal 
in the equatorial regions, while being far more signihcant in the polar regions. In the 
equatorial regions, the A||S for each of the Parker HMF, SBM and JKM are identical at 
small r and then start diverging from one another at about 1 AU. This divergence is not 
signihcant however, and both the A||S for the SBM and JKM end up being only marginally 
below the Parker HMF, the SBM being the lowest. In the polar regions, the A||S for each 
of the Parker HMF, SBM and JKM are again seen to start out together at small r and 
then to diverge at greater r. As should be the case, the A||S for the SBM and JKM are 
observed to start diverging from the Parker HMF at more or less the same radial distances 
at which the HMF prohles in Figure started diverging from one another. In the polar 
regions the A||S for both the SBM and JKM then keep signihcantly lower than that for 
the Parker HMF, the JKM being the lowest up to r = tts where the SBM drops below it 
due to the failure of the JKM to decrease over the TS. 


The drift coefficient also has a 1/B dependence according to Eq. (13), and a modih- 


cation to the HMF will therefore ahect the particle drift velocity, so that in each of the 
cases for the JKM and SBM a set of expressions analogous to Eq. (16) can be obtained. 


For the modihed JKM, i.e. using Eq. (30), the drift velocities of Eq. (16) are replaced by 


Vdr = --ydo I 1 + (jr^ ^ d' cot 6^ 
Vd9 = -ydo I 2 + j T 

Vd4> = Udo ( 4f^cot6'^ ^2-h 


(33) 
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with 


■ydo — 2 - 


K,. 


rB\ 


'l+f^l +4-^ 

Vq, 


-3 • 


(34) 


To obtain the analogous set of equations in the case of an HMF modified according to 
the SBM, we write 

^ ^ - 6'r, (35) 


where the SBM modification term 5' is dehned as 

O.O2f"sw(&,0) 
b f"sw(r,0)’ 

since we have set B-^ih)/Bj^ih) = —0.02 for our choice of 6; hnally, substituting 

'y = rS — 4/, 


(36) 


(37) 


Eq. (16) is replaced by 


Vdr = (7^ + 7^^ + 7 - d') cot 0 

Vde = -2vdo (2 + 7 ^) 7 
Vdy, = -2ndo7^ cot 6 , 


with 


^'do = 


Ka 


(38) 


(39) 


tSa/I + 7^ 

Figure shows the drift velocity streamlines in the case of the SBM. Compared to 
Figure it is clear that the inward drift from the poles is not as effective when employing 
the SBM as it was when using the Parker HMF: the lines in the case of the SBM fan out 
significantly more than in the case of the Parker HMF. The value of B^^b)/Bji{b) = —0.1, 
and has once again been set to a value low enough so that these differences can be clearly 
illustrated. 


4.4 The Effect on Modulation 

Very effective use can be made of the properties of the SDE-based model to qualitatively 
illustrate the effect that these modihcations to the HMF can have. Figure presents a 
meridional cut of the heliosphere and shows trajectories for pseudo-particles (1.5 GeV 
protons) entering in the polar regions during an H > 0 polarity cycle, and propagating 
down towards the Sun, which is located at (Y, Z) = (0, 0); the solid black line shows the 
warped HCS with tilt angle a = 10°, for illustrative purposes. The position of the TS, 
tts = 88 AU, is indicated by the vertical red lines. The red, blue and green trajectories 
respectively shows the case for the Parker HMF, JKM and SBM, and the differences are 
signihcant: in the case of the Parker HMF, the drift is obviously more effective, followed by 
consecutively reduced drift effects, i.e. reduced in terms of the ratio to diffusive processes, 
in the case of the JKM and SBM. This is not only evident from the longer trajectories in 
the latter two cases, but also by the reduced occurrence of the long straight and jump- 
like segments in the trajectories, which indicate cases where Vd reaches unphysically large 
speeds. From this hgure, therefore, it is evident that both the JKM and SBM are effective 
in reducing drift effects over the poles and are thus successful in their original purpose. 
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Next, we look at the effect that these HMF modifications have on CR proton spectra; 
in this way a quantitative measure of the effect on the modulation is obtained. The model 
is used to reproduce the December 2009 CR proton spectrum, which was observed by the 
PAMELA experiment and averaged over 1 month (e.g. Adriani et ai, 2013; Potgieter 
et ai, 2014). This is shown in Figure where the observations from PAMELA are 
indicated by the circles and the model reproduction, assuming the SBM, is indicated by 
the green line; the LIS at 120 AU is represented by the black line. The values of all the 
relevant quantities for this model reproduction are indicated in the last column of Table 1. 
To obtain a comparison with the SBM, these same values are used while assuming the 
Parker HMF (red line) and JKM (blue line) respectively; the difference between these 
three HMF profiles is clearly depicted: the JKM and SBM are consecutively lower than 
the Parker HMF; these spectra peak in the vicinity of 200 MeV and, at this energy, the 
SBM intensities are reduced relative to the LIS by about 86%; the JKM intensities are 
reduced by about 82% while, in the case of the Parker HMF, intensities show a decrease 
of only about 71%. The fact that the spectrum resulting from the SBM is lower than 
that resulting from the JKM once again illustrates the fact that, although the effects of 
the JKM may be locally greater than that of the SBM, the overall effect of the SBM is 
larger than that of the JKM when integrated over the whole of modulation space. 

Further motivation for the applicability of the SBM comes from applying this mod¬ 
ification to model the peculiar solar minimum stretching from 2006 to 2009. Monthly 
averaged proton spectra measured by PAMELA at Earth, during November 2006 and 
December of 2007 and 2008 are now also reproduced using our model. This is depicted 
in Figure IR, the observations from PAMELA once again indicated by the circles, while 
the model reproductions for 2006, 2007, 2008, and 2009 are indicated by the red, orange, 
green, and purple lines respectively; the LIS at 120 AU is represented by the black line. 
The parameters to obtain these results are summarised in Table 1 and it is clear that: 1) 
the diffusion was required to continually increase from 2006 to 2009, with the value of K||,o 
increasing from 16.0 in 2006 to 17.1, 17.5, and 19.4 in 2007, 2008, and 2009 respectively; 
2) the spectra became progressively softer towards 2009, so that the rigidity dependence 
of the diffusion coefficient below ~ 3 GeV was required to change from in 2006 to 
pO-81^ po.78^ po.73 2007, 2008, and 2009. These type of changes are in agreement with 

those found by various other authors, including Bazilevskaya et al. (2012), Ndiitwani et 
al. (2013), Potgieter et al. (2014), and Pacini and Usoskin (2015). Except for providing 
evidence for the applicability of the SBM in particular, this qualitative correspondence in 
modulation results over the entire solar minimum period of 2006 to 2009 also illustrates 
the already well established credibility of the SDE-based model. 

Figure 11 illustrates the results of Table 1 by depicting the mean free paths for each 
of 2006, 2007, 2008, and 2009. The grey lines show the mean free paths obtained by 
Potgieter et al. (2014) and, when compared to the results obtained in this study, the 
qualitative agreement is clearly illustrated. Quantitative differences must exist between 
these two sets of results for three reasons. Firstly, these two studies did not select the same 
HMF profile, i.e. this study employs the SBM, while the study of Potgieter et al. (2014) 
utilised the JKM; referring to Figure this will clearly have a quantitative influence on 
the results. Secondly, our model uses an updated LIS that was adapted according to 
observations by Voyager 1 after its crossing of the HP, as explained earlier. Thirdly, there 
were some differences in the treatment of current sheet drift, the details of which are not 
relevant to the topic of this paper. 
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Table 1: Summary of the parameters used in this study to reproduce the 2006 to 2009 
PAMELA proton spectra, using the SBM; see Eqs. 0, (|^ and (31). 


Parameter 

2006 

2007 

2008 

2009 

a [deg] 

15.7 

14.0 

14.3 

10.0 

Be [nT] 

5.05 

4.50 

4.25 

3.94 

^TS 

88.0 

86.0 

84.0 

80.0 

Kll^o [6 X 10^° cm^.s“^] 

16.0 

17.1 

17.5 

19.4 

ai 

0.85 

0.81 

0.78 

0.73 

Po [GV] 

ol 

1/^ 

1/v^ 

1/yi) 

a [GV] 

4.0 

4.0 

4.0 

4.2 


5 Conclusions 

An investigation into the effects of the SBM and JKM modihcations were undertaken, 
making use of an SDE-based numerical modulation model. These two modihed HMF 
prohles were compared to that of the unmodihed Parker HMF and it was illustrated that 
both these modihcations change the radial dependence of the HMF in the polar regions, 
such as to bring about a larger HMF magnitude as function of radial distance. This 
larger HMF magnitude in the polar regions then reduces both the diffusion and drift 
coefficients in these regions. It was noted that, although the effects of the JKM may be 
locally larger than those of the SBM, the effects of the SBM are exercised over a greater 
part of modulation space and hence its overall effect on CR intensities is greater than 
that of the JKM. The eventual effect on modulation of these two modihcations were also 
illustrated and it was concluded that both are effective in reducing drift effects over the 
poles; this was illustrated by plotting drift velocity streamlines, as well as utilising the 
unique abilities of the SDE-based numerical model, which enables the construction of 
pseudo-particle trajectories. It was then shown that the SBM can be used to reproduce 
the Galactic proton spectra observed by the PAMELA experiment during the peculiar 
solar minimum period of 2006 to 2009, and that the results thus obtained correspond 
qualitatively to that found by various other authors who have previously investigated this 
solar minimum. This agreement was presented as veritable evidence for the applicability 
of the SBM in CR modulation studies; it was also argued that the SBM is grounded in 
observational evidence and that, for these reasons, it should be preferred as modihcation 
to Parker type HMFs. 
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Figure 1: Magnetic field lines for the case of the Parker HMF, originating from consecutive 
azimuthal angles and a latitude of 9 = 45°; shown in Cartesian coordinates and over the Erst 
10 AU from the Sun, which is located at (X,V,Z) = (0,0,0); the colour scale corresponds to 
r‘^\B\. The spiral structure of this held is clearly depicted. 



Figure 2: Gradient and curvature drift streamlines in the case of the Parker HMF, shown in 
Cartesian coordinates. From the narrow, cone-like conEguration of these lines it is evident that 
drift through the polar regions of the heliosphere proceeds effectively. The Sun is again located 
at {X,Y,Z) = (0,0,0). The colour scale corresponds to the magnitude of the gradient and 
curvature drift velocity. 
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Figure 3: Magnetic field lines resulting from employing the SBM with B^{h )/= —0.5. 
Compared to Figure^ these lines are evidently closer together than in the case of the unmodified 
Parker HMF, in accordance with the larger HMF magnitude of the modified held. The held 
lines again start from successive azimuthal angles, latitude 9 = 45°, and is shown over the hrst 
10 AU from the Sun; the colour scale again corresponds to r‘^\B\. 
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Figure 4: Magnetic field magnitude radial profiles in the equatorial regions (top panel) and in 
the polar regions (bottom panel). The Parker HMF is indicated by the solid red line while the 
SBM and JKM are represented by the dashed green and blue lines respectively. 
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Figure 5: Magnetic field latitudinal profiles at four selected radial distances r = 1 AU (first 
panel), 10 AU (second panel), 50 AU (third panel), and 100 AU (fourth panel). The Parker 
HMF is indicated by the solid red line while the SBM and JKM are represented by the dashed 
green and blue lines respectively. 
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Figure 6: Radial profiles for the parallel mean free paths A|| of a 1 GV proton for each of the 
Parker HMF (red lines), the SBM (green lines), and the JKM (blue lines). The upper three lines 
depict the A||S in the polar regions, while the lower three lines do so in the equatorial regions. 



Figure 7: Gradient and curvature drift streamlines, upon employing the SBM with 
B'j'{b)/B^{b) = —0.1. Compared to Figure^ the cone-like configuration of these lines are 
far less narrow, i.e. it is significantly more fanned out. This indicates less effective inward drift 
than in the case of the Parker HMF. Once again, the colour scale corresponds to the magnitude 
of the gradient and curvature drift velocity. 
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Figure 8: Pseudo-particle trajectories reflected in the meridional plane of the heliosphere for 
1.5 GeV protons undergoing drift and diffusion during an A > 0 drift cycle where a = 10°. 
The black line depicts the wavy HCS, while the red, green and blue lines respectively show the 
pseudo-particle trajectories through the polar regions in each of the cases for the Parker HMF, 
the SBM, and the JKM. The vertical red lines are drawn at the position of the TS, rxs = 88 
AU. The Sun is located at {Y, Z) = (0, 0) with Y the equatorial direction. 
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Figure 9: The SBM (green line) is applied to reproduce the 2009 Galactic proton spectrum as 
measured hy PAMELA (circles) at Earth and published in e.g. Potgieter et al. (2014). The 
blue and red lines respectively show the effects of the JKM and Parker HMF, using the same 
parameters as for the SBM. The black line represents the LIS at r = 120 ALf. 
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Figure 10: The circles represent the 2006 to 2009 Galactic proton spectra from observations, at 
Earth, by the PAMELA mission, as published by e.g. Potgieter et al. (2014); these spectra, from 
the lowest to the highest, correspond to November of 2006 and December of 2007, 2008, and 2009 
respectively. The SDE-based model of this study is applied, utilising the SBM, to reproduce 
these observed spectra; the reproductions for 2006, 2007, 2008, and 2009 are indicated by the 
red, orange, green and purple lines respectively, and quite clearly become higher as well as softer 
toward 2009, with respect to the LIS (black line) at r = 120 AU. 
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Figure 11: The parallel mean free paths A||S for 2006, 2007, 2008, and 2009 obtained from the 
SDE-based modeling shown in Figure are represented by the red, orange, green, and purple 
lines respectively. For comparison, the A||S that resulted from the modeling of Potgieter et al. 
(2014) are included and represented by the grey lines, so that the qualitative correspondence 
between these two sets of results are clearly depicted. 
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